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Abstract 

We study the phenomenon of negative energy densities in quantum field theories with self¬ 
interaction. Specifically, we consider a class of integrable models (including the sinh-Gordon model) 
in which we investigate the expectation value of the energy density in one-particle states. In this 
situation, we classify the possible form of the stress-energy tensor from first principles. We show 
that one-particle states with negative energy density generically exist in non-free situations, and 
we establish lower bounds for the energy density (quantum energy inequalities). Demanding that 
these inequalities hold reduces the ambiguity in the stress-energy tensor, in some situations fixing 
it uniquely. Numerical results for the lowest spectral value of the energy density allow us to 
demonstrate how negative energy densities depend on the coupling constant and on other model 
parameters. 


1 Introduction 


The energy density is one of the fundamental observables in classical as well as quantum field theories. 
It has a special significance in field theories on curved backgrounds, since it enters Einstein’s field 
equation as a source term, and is therefore linked to the structure of space-time. But also on fiat 
Minkowski space, as well as in low-dimensional conformal field theories, it plays an important role. 

In the transition from classical to quantum field theories, some of the distinctive properties of the 
energy density are lost. In particular, the classical energy density is positive at every point, which in 
General Relativity implies certain stability results, such as the absence of wormholes |FR96j . This is 
however not the case in quantum field theory, not even on fiat spacetime: While the global energy 
operator H is still nonnegative, the energy density can have arbitrarily negative expectation values 
|EG.I65j . However, a remnant of positivity is still expected to hold. When one considers local averages 
of the energy density, = f dt g^{t) 0) for some fixed smooth real-valued function g, then 

certain lower bounds - quantum energy inequalities (QEIs) - should be satisfied. In the simplest case, 
one finds for any given averaging function g a constant Cg > 0 such that 

{ip,T^^{g‘^)(p) >-CgWiff ( 1 . 1 ) 


for all (suitably regular) vector states (p of the system; this is a so-called state-independent QEI. In 
general, (1.1) may need to be replaced with a somewhat weaker version (a state-dependent QEI) where 
the right-hand side can include a slight dependence on the total energy of the state (p. 

This raises the question under which conditions the QEI (1.1) can be shown to hold rigorously. 
The inequality has in fact been established for various linear quantum fields, on fiat as well as curved 
spacetime, and in conformal QFTs (e.g., |Fla97(lPF98tIFewOnlIFVn2lIFHDb] : see |Fewl2| for a review). 
However, dropping the restriction to linear fields, that is, allowing for self-interacting quantum field 
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theories, few results are available. This is not least due to the limited availability of rigorously 


constructed quantum field theoretical models; see however our recent proof of (1.1) in the massive 


Ising model |BCF13] . In a model-independent setting, one can establish state-dependent inequalities 
for certain “classically positive” expressions |BFn9| . based on a model-independent version of the 
operator product expansion |Bos05j , but the relation of these expressions to the energy density remains 
unknown. 


In the present paper, we will investigate the inequality (1.1) in a specific class of self-interacting 


models on 1-|-1 dimensional Minkowski space, so-called quantum integrable models, which have recently 
become amenable to a rigorous construction. Specifically, we consider integrable models with one 
species of massive scalar boson and without bound states. 

An a priori question is what form the stress-energy tensor takes in these models. There 
is a straightforward answer in models derived from a classical Lagrangian, such as the sinh-Gordon 
model, where a candidate for the operator can be computed |FMS931 IKM931 [MS94] . However, we also 
consider theories where no associated Lagrangian is known; and more generally, we aim at an intrinsic 
characterization of the quantum theory without referring to a “quantization process”. In fact, there 
will usually be more than one local field that is compatible with generic requirements on T“^, such as 
covariance, the continuity equation, and its relation to the global Hamiltonian. 

Given one can ask whether the QEI (ll.ll) holds for this stress-energy tensor, or rather for 


which choice of stress-energy tensor. In fact, QEIs may hold for some choices of but not for 
others, as the nonminimally coupled free field on Minkowski space shows |EO08) . Ideally, one would 
hope that requiring a QEI fixes the stress-energy tensor uniquely. 

In the present article, we will consider the above questions in integrable models, but with one 
important restriction: We will consider the energy density at one-particle level only. That is, we will 


ask whether the inequality (1.1) holds for all (sufficiently regular) one-particle states (p. 

This case might seem uninteresting at first: One might argue that in integrable models, where 
the particle number is a conserved quantity, the effect of interaction between particles is absent at 
one-particle level. But this is not the case, as already the massive Ising model shows |BGF13j : in this 
model of scalar bosons, one-particle states with negative energy density exist, whereas in a model of 
free bosons, the energy density is positive at one-particle level. We will demonstrate in this paper that 
self-interaction in our class of models leads to negative energy density in one-particle states, and that 
this effect increases with the strength of the interaction. 

Our approach is as follows. Having recalled the necessary details of the integrable models con¬ 
sidered (Sec. [^, we ask what form the stress-energy tensor can take at one-particle level. This will 
lead to a full characterization of the integral operators involved, with the energy density fixed up to 
a certain polynomial expression (Sec. [^. 

In Sec. 1^ we will show that under very generic assumptions, states of negative energy density 
exist, and that for certain choices of the stress-energy tensor, the energy density becomes so negative 
that the QEI ( |1.1[ ) cannot hold even in one-particle states. For other choices of the energy density 
operator, we demonstrate in Sec. [^that the QEI does hold. In consequence (Sec. [^, we find in the 
massive Ising model that one-particle QEIs hold for exactly one choice of energy density, whereas in 
other models (including the sinh-Gordon model), the choice of energy density is at least very much 
restricted by a QEI. 

All this is based on rigorous estimates for the expectation values of However, the best 

possible constant Cg in (1.1) - in other words, the lowest spectral value of restricted to one- 


particle matrix elements - can only be obtained by numerical approximation. We discuss the results 
of an approximation scheme in Sec. thus demonstrating how the effect of negative energy density 
varies with the coupling constant and with the form of the scattering function. The program code 
used for this purpose is supplied with the article |Godj . We end with a brief outlook in Sec. 

2 Integrable models 


Eor our investigation, we will use a specific class of quantum field theoretical models on 1-|-1 dimen¬ 
sional spacetime, a simple case of so-called integrable models of quantum field theory. These models 
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model 

5(C) 

FminiC) 

parameters 

free field 

1 

1 


Ising 

-1 

—i sinh 1 


sinh-Gordon 

sinhC — isinB7r/2 
sinhC + isini?7r/2 

exp Jb(C) 

0 < B < 2 

(generalized 

sinh C, — i sin Bji: j‘1 

exp Ei=i«^B,(C) 

Bj G (0, 2) -|- zM, 
either real or in 

sinh-Gordon) 

^ sinh C + * sin Bj7r/2 

(—i sinh i)2L«./2J 

(generalized Ising) 

sinh (f — i sin Bj7r/2 

exp Jb, iO 

complex-conjugate 

sinhC + isinBjTT/2 

j = l ^ ■! 

{—i sinh k) 2 rn/ 2 ]-i 

pairs 


Table 1: Examples of scattering functions S and corresponding minimal solutions -Fmin 
for the definition of the function Jb- 


See Eq. 


(2.5) 


describe a single species of scalar massive Bosons with nontrivial scattering. The scattering matrix is 
factorizing: When two particles with rapidities 6 and r] scatter, they exchange a phase factor S{9 — rj)] 
and multi-particle scattering processes can be described by a sequence of two-particle scattering pro¬ 
cesses. 

There are several approaches to constructing such integrable quantum field theories. Convention¬ 
ally, one starts from a classical Lagrangian, derives the two-particle scattering function S from there, 
and then constructs local operators (quantum fields) by their matrix elements in asymptotic scatter¬ 
ing states; this is the form factor programme |Smi92llBF!ai6] . A more recent, alternative approach 
|SW00| ILecOHI starts from the function S as its input, then constructs quantum field localized in 
spacelike wedges (rather than at spacetime points), and uses these to abstractly obtain observables 
localized in bounded regions. 

We will largely follow the second mentioned approach here; in particular, we set out from a function 
S rather than from a classical Lagrangian. In all what follows, we will assume that a scattering function 
S is given, which we take to be a meromorphic function on C which fulfils the symmetry properties 

5(-C) = 5(C)-' = 5(C + ivr) = W)- (2.1) 

A range of examples for such functions 5 can easily be given, in particular because the properties 
are preserved under taking products of functions; see Table This includes the sinh-Gordon model, 
depending on a coupling parameter B, which is normally constructed from a Lagrangian |EMS93| : but 
for other examples (e.g., the generalized sinh-Gordon models mentioned in Table[^, no corresponding 
Lagrangian is known. 

We will not enter details of the construction of the associated quantum field theory based on 5 
here, but will recall only the general concepts as far as relevant to the present analysis. The single 
particle space of the theory is given by Tii = L‘^(R,d9), where the variable 9 of the wave function is 
rapidity, linked to particle two-momentum p by p{9) = /r(cosh 9, sinh 9 ); here /r > 0 is the particle mass. 
On T-Li, the usual representation of the Poincare group acts. One then constructs an “5-symmetric” 
Eock space T-L over Tii, on which “interacting” annihilation and creation operators z{9) and z^{9) act; 
instead of the GCR, they fulfil the Zamolodchikov-Faddeev relations [LecOSj . depending on 5. The 
quantum field theory is constructed on this Eock space. If A is an operator of the theory (of a certain 
regularity class, including smeared Wightman fields), localized in a bounded spacetime region, then 
it can be written in a series expansion |BG13l IBGlSj 

A= ^ / ^^^^Fl^l^{0 + iO,r] + i7v-iO)z\9i)---z\9„i)z{r]i)---z{r]n), (2.2) 

m,n=0 

where Fm+n are meromorphic functions with certain analyticity, symmetry and growth properties 
(which we will recall where we need them). Examples of such local observables A would include 
smeared versions of the energy density, supposing they fall into the regularity class mentioned. 
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In the construction of functions that fulhl these properties, an important ingredient is the so- 
called minimal solution -Fmin of the model [KW78) . We consider it here with the following conventions. 

Definition 2.1. Given a scattering function S, a minimal solution is a meromorphic function i^min 
on C which has the following properties. 

(a) Fmin has neither poles nor zeros in the strip 0 < Im(^ < vr, except for a first-order zero at ( = 0 
in the case that 5(0) = —1, 

(b) FmmiiTT + C) = “0/ 

(c) ih[iin(~C) = 5( —^)i^min(C)) 

(d) Fmin(i7r) = 1; 

(e) There are constants a,b > 0 such that | log|Fmin(C)| | < a|ReC| +b i/|ReC| > 1, 0 < Imi) < vr. 


Note that properties (a) and (e) automatically hold analogously for the strip M -|- z[7r, 27r] by 
property The hrst-order zero at C = 0 (and analogously C = ^vr) must necessarily occur in the 

case 5(0) = —1 due to ( 


The properties 


actually fix Tmin uniquely if it exists, so that we can speak of the minimal 


solution. We prove this in our context; cf. [KW78t p. 459]. 



G(C + 27ri) = G(C) = G(-C). 


(2.3) 


We can therefore find an entire function P such that G{C) = P(cosh((): We observe that cosh( •) 
is bijective from the region (0, oo) ± i(0,7r) to the upper, respectively lower, half-plane, and use the 


properties (2.3) to accommodate the branch cuts of the inverse hyperbolic function. From property 


|(d)[ this function fulhls the estimate 

log|P( 2 ;)| < a'|Rearcoshz| -|-6' 


(2.4) 


with certain constants a',b' > 0 and for large Rearcoshz. Since Rearcosh^ grows like log2|2;| for 
large \z\, this means that P is polynomially bounded at infinity, and hence a polynomial. But due to 
property (a), P has no zeros, and is therefore constant. Now from [(d)| P{z) = P(—1) = 1. Hence 

□ 


7?^ — 

'^min '^min 


Let us note a s imple consequence: One checks that together with Pminj also Cjz]’ -^min(—C) fulfils 

Together with (b), this shows that 


properties (a) - (e) Thus the lemma yields Pr 
Pnin is symmetric and real-valued on the line I 


i( C) — Pnin(C)' _ 

-|-i7r. We will use this fact frequently in the following. 
In this article, we will always assume that a minimal solution Pmin exists. In fact, for the examples 
we mentioned, they are listed in Table For the sinh-Gordon and related models, this involves the 
integral expression 


Jb{C + *7r) := 8 


dx sinh ^ sinh sinh § 


2 ■ 2 


sinh^ X 


sm 


— for R G (0, 2) -|- il 

ZTT 


(2.5) 


which is known from |FMS93j (but note that we use a different normalization for Pmin)- Due to the 
Riemann-Lebesgue lemma, Jb{0 converges to a constant as |ReC| —^ oo, with 0 < ImC < 27r fixed. 

For the generalized sinh-Gordon and Ising models in Table Pmin can essentially be obtained as 
a product of the minimal solutions of the corresponding sinh-Gordon or Ising factors, since properties 


(b) - (e) in Def. 2.1 are again preserved under products. However, in order to satisfy property (a), any 
possible double zeros of the product function at C = 0 need to be cancelled by dividing by appropriate 
powers of (—isinh(C/2))^. 
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3 Energy density at one-particle level 


The first question we want to consider is which form the energy density operator can take in our 


models. More specifically, we ask what the functions in the expansion (2.2) can be if 


yl = r“^(/)= / dt/(t)r“^(t,0) 


(3.1) 


is a component of the stress-energy tensor smeared with a real-valued test function in time direction. 
The answer may appear obvious in models such as the free field or the sinh-Gordon model, where the 
energy density is linked to the classical Lagrangian and well studied. However, we aim at an intrinsic 
characterization of the energy density within the quantum theory, and therefore we are looking for the 
most general form of the stress-energy tensor compatible with generic assumptions on this operator, 
which will be detailed below. 

As announced in the introduction, we will consider only one-particle states of the theory, and 
evaluate the stress-energy tensor only in these. More precisely, we will consider the stress-energy 
tensor only in matrix elements of the form 


((^, r"^(/) V') with (^,'0 G T>(M) C ?^ 1 . 


(3.2) 


(The restriction of the quadratic form to smooth functions of compact support, i.e., V’ £ ^(l^)j is 
perhaps too cautious - the form can easily be extended to non-smooth and to sufficiently rapidly de¬ 
caying wave functions, and we will in fact use piecewise continuous functions in the numeric evaluation 
in Sec. but for the moment we restrict to 2?(M) for simplicity.) 

Writing in expanded form as in (2.2), we see that only the coefficients To and F 2 contribute 

to this one-particle matrix element. Since the coefficient To[A] equals the vacuum expectation value 
of the operator, we can assume without loss that To[r"^(/)] = 0; in fact, this is necessary for the 
energy density, since it would otherwise not integrate to the Hamiltonian H. This leaves us with 


{(f, = j de dr] (p{e)F2[T°‘^ {f)]{e + i0,ri + i7T - iO)V'(??). (3.3) 

Since we expect to be a translation-covariant operator-valued distribution in /, we will assume 

that 

F 2 [T°‘^{f)]{9, T] + irr) = F^^{6, r])f[]i cosh. 0 — ]i coshrj) (3.4) 

with a function independent of /, where we take the Fourier transform with the convention 
f{p) = f dt f{t)F^P. 

The task is therefore to determine the possible form of starting from physical properties. We 
will list these assumptions one by one and explain their motivation, but we skip details of how they 
are derived; that is, we will take these assumptions as axioms in our context. 

The first set of conditions follows from the general properties of the expansion coefficients 
of a local operator A in integrable models, as derived in |BC15] . 


(Tl) are meromorphic functions on C^, analytic in a neighbourhood of the region —vr < lm(0 — 

7 ]) < TT. 

This is due to the general analyticity properties of the coefficients F^ |BC151 property (FDl)], 
together with the absence of “kinematic poles” in the specific case of the coefficient F 2 (from 
property (FD4) there). 

(T2) They have the symmetry properties 

F^^{e,ri) = S{e -r])F^^{ri -\- ITT, 6 — in) = F°^^{t]-\- in, 6 + in). (3.5) 


This is a rewritten form of the properties of “S'-symmetry” and “S'-periodicity” - properties 
(FD2) and (FD3) in |BG15j . 
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(T3) Hermiticity of the observable is expressed as 


rj) = F°‘^{rj, 9) for all 0, r/ G M. 


(3.6) 


(T4) We demand that there exist constants k,i > 0 such that 

|F“^(0,r/)| < £(coshRe6l)^(coshRer/)^ whenever — vr < lm(0 — r]) < it. (3.7) 


This is motivated by the bounds for Fk [^] discussed in [BC15] ; the condition guarantees that the 
smeared version T^^[f) for every Schwartz function / will fulfil, at one-particle-level, polynomial 
high-energy bounds in the form of |BC15l property (FD6)]. 


Further conditions are derived from properties that one expects specifically of the stress-energy tensor 

'j'afS 

(T5) The stress-energy tensor is a symmetric tensor, which rewrites in our context as 

F"^(0,r/) =F^“(0,r/). (3.8) 


(T6) It is covariant under Lorentz transformations, which means in our terms, cf. [BC131 Prop. 3.9], 

F^P{0 _ A, 7/ - A) = r/), (3.9) 

where A is the boost matrix with rapidity parameter A. 

(T7) It is invariant under spacetime reflections, which by |BC15L Thm. 5.4] translates to 

F“^(0 + i7r,77 + i7r) =F“^(0,r/). (3.10) 


(T8) It fulhls the continuity equation [daT^^ = 0), which reads in our terms, 

{PaiO) - Pa{v))F'"^i0, v) = 0, (3.11) 

where p{9) = ;u(cosh0, sinh0). 

(T9) The (0,0)-component of the tensor integrates to the Hamiltonian, f dxT^^(t,x) = H, which 
translates to ^ 

F^^{d,e) = ^coshH. (3.12) 

27r 

Taking these properties as our starting point, we ask what functions F°‘^ are compatible with 
them. The answer is given in the following proposition. 

Proposition 3.1. Functions fulfil the properties 
polynomial P with P{1) = 1 such that 

F“^(0, 7?) = F“£(0, 7y)P(cosh(0 - 7/))F^i„(0 - 7? + ivr), (3.13) 


(T1)- (T9) if, and only if, there exists a real 


where 




Ijf f cosh^(^^) ^ sinh(6*-I-7/)A 

27 r sinh(0-b 77) sinh^(^^) / 


(3.14) 


Note that 
the free Bose field. 


is the well-known one-particle expression of the “canonical” stress-energy tensor of 
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Proof. It is straightforward to check that as given in (3.13) fulfils all conditions (T1)-(T9) 
knowing that F^^^ fulfils them in the case S = 1. 


Thus, let F°‘^ fulfil conditions (Tl) - (T9) 
with |(T5)| to obtain 


We first use 


(T 8 ) with (3 = 0 and with (3 = 1, along 


F^^{9, rf) = tanh 


2 ^ ^ pOO 


F^\e,r^). 


(3.15) 


Now we consider the functions G“^(C) := F°^^{((/2, —((/2), which are meromorphic by (Tl) Note that 
|(T 8 )| with (3 = 0 implies 

(po(C) - Po(-c)) G°0(2C) + (pi(C) - Pi(-C)) G10(2C) = 0 (3.16) 


=0 


=—2/i sinh Q 


and hence = 0. Then = 0 by (T5) Also, (3.15) leads to = 0, and the only nonzero 

(3.17) 

(3.18) 


component of G is G^^. From (T2) and (T7) one concludes 

G“(-C) = GOO(C), G“(C - ivr) = S(C) G“(C + ivr). 

Now set 

Ott „„ 

i(C + ivr). 


QiO := 5g»»(C)/F. 


This Q is analytic in a neighbourhood of the strip M + i[—vr, vr]. (Note that in the case S'(O) = —1, 
the zeros of the denominator at C = TivT are cancelled by corresponding zeros of the numerator which 
exist due to (3.17).) The symmetry relations (3.17) and Def. ^ imply that 

Q{-C) = Q{C), Q(C + * 7 r) = Q(C-ivr). (3.19) 


Arguing as in the proof of Lemma 2 . 2 , we can therefore find an entire function P such that Q{C) = 
P(cosh(^). From property |(T4)| and Def. 2.]|^, this P fulfils the estimates 

27r |G°‘^(arcosh 2 :)| ^ ^, (cosh 7 Re arcosh z) 


\Piz)\ < 


<^'( 2 |z| + l) 


k—k' 


(3.20) 


|Tmin(arcosh(z) + Avr)! (coshRe arcosh 2 ;)*^' 

with some con stants £' ,k' . That is, P is a polynomially bounded entire function, and hence a poly¬ 


nomial. From 


(T3), we can conclude that P{z) = P{z) for all 2 , thus the coefficients of P are real. 


Also, (T9) implies G*^®(0) = /r^/27r and hence P(l) = 1. Thus P has the properties claimed in the 
proposition. Combining our results, we have shown that 

^ (J 0) P(coshC)Pmin(C + ivr). 


(3.21) 


This means that (3.13) holds in the case 6 + r] = 0. But then it holds for any value ol 9 + rj, since 
both sides of (3.13) fulfil the covariance condition |(T 6 )[ □ 

In the following, when we speak of the stress-energy tensor of a model, we will always refer to one 
of the form (3.4), with as in Proposition 3.1 We will often abbreviate 


FpiC) ■= P(cOshC)Pmin(C + ^Tt), 

noting that Pp enjoys most of the defining properties of Pmin (namely, Def. 2.1 'b) 


(3.22) 
with shifted 

argument, but not (a)). Also, Pp is symmetric and real-valued on the real line. The expectation value 
of the energy density now becomes 

,2 _ 

(3.23) 






‘ _ Q I ^ 

d9 dri (p{9)(p{r]) cosh^ —^ Fp{9 — rf) /{fj. cosh 9 — //coshr/). 

Proposition |3.1| shows that on the one-particle level, we recover the well-known “canonical form” 
of the energy density of the free field and of the sinh-Gordon model |FMS93j , and the considered 
for the massive Ising model in [BCF13j . up to a possible polynomial factor P in cosh(0 — 77 ). We 
emphasize that, while one might expect P = 1, all our assumptions so far are perfectly compatible 
with a more general polynomial P. However, we will see later (in Sec. that the choice of P is 
restricted, in some cases uniquely to P = 1, if we demand that quantum energy inequalities hold. 
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4 States with negative energy density 


As the next question about properties of the energy density, we will ask whether single-particle states 
with negative energy density exist at all; more specifically, whether can be negative if 

(p G 2?(M) and g' is a real-valued Schwartz function, i.e., g G 5r(M). The example of the free held with 
canonical energy density {S = 1, Fp = 1) shows that this is not guaranteed: In this specihc case, 
T^^[g‘^) is known to be positive between one-particle states. However, as we shall see in a moment, 
the introduction of interaction quite generically leads to negative energy densities. 

We will exhibit these negative energy densities by explicitly constructing corresponding states p. 
In preparation, we hx a nonnegative, smooth, even function y with support in [—1,1]. For q > 1, 
p > 0 , we set Xg,p(6') := p~'^^‘^\\x\\q^x{0/p), so that Xq,p has support in [—p, p] and is normalized with 
respect to the norm || • ||q. 

Proposition 4.1. Suppose that there is 6p £ U. such that \Fp{6p)\ > 1. Then there exist g G 5r(M) 
and p G P(M) such that {p,T^'^{g‘^)p) < 0. 


Proof. We will show that, with suitable choice of p, one has 

0 + r] 


0 > dOdp p{9)p{r]) cosh^ “ ' ' Fp{9 — p) =: X. 


(4.1) 


{X is the expectation value of r^^(O) up to a factor.) Rewriting (3.23) as 

,2 


{p,T^\g^)p) = ^j dtg\t) I d 0 dr? 99 ( 0 )¥p(r?)cosh 2 ^Fp( 0 -r?)e**^(“«‘'^-“^^^), (4.2) 

and noticing that the inner integral expression is real, continuous in t, and gives A < 0 at t = 0 , it is 


then clear that we can choose g so that (4.2) becomes negative. 
To achieve (4.1), we will choose the wave function p as 

2 

= -Ij), 

i=i 


(4.3) 


where fdj G C, Xj ^ p > 0 will be specified later; the quantity X then depends on these 

parameters. To show that A < 0 for some /? > 0, it suffices to show that A converges to a negative 
limit as p 0. Noting that xi,p{0) ^(0) this limit, one obtains from (4.1) that 


lim A = 3*M3, where 

p\0 


Mjk 


= cosh^ - Ik)- 


(4.4) 


This expression is negative for suitable /3 = (/ 3 i,/ 32 ) if the determinant of the matrix M is negative. 
Setting 7 i := 7 -|- 0 p/ 2 , 72 1=7 — 0p/2, with 7 G M still to be chosen, one computes 

detM = cosh^( 7 -|-0p/2) cosh ^(7 — 0p/2) — cosh^ 7 Fp( 0 p). (4.5) 

Since cosh (7 -|- 0p/2) cosh (7 — 0pf2)l cosh^ 7 —)■ 1 as 7 —)• 00 , and since Fp{0p) > 1 by assumption, 
det M does indeed become negative for sufficiently large 7 , which concludes the proof. □ 


In particular, the conditions of Proposition 4.1 are met in the sinh-Gordon and the Ising models 
for any choice of P, as well as in the free model if P ^ 1. Thus single-particle states with negative 
energy density exist in generic situations. 

Under stricter assumptions on the function Fp, we can in fact show a significantly stronger result: 
If Fp grows stronger then a certain rate, then the negative expectation values of T^^{g‘^) become so 
large that quantum energy inequalities cannot hold. 

Proposition 4.2. Suppose there exist 00 > 0 and c> ^ such that 

^0>0o: Fp{0) >ccosh0. (4.6) 

Let g G 5r(M), g ^ 0. Then, there exists a sequence {pj)j^-^^ in 'P(M), \\pj\\ = 1, such that 


{9 )Tj) -00 as j 00 . 


(4.7) 
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Proof. We set 


(Pj{9) := 1 X2,pj {0 - j) + Pj,2 X2,pj (9 + j), (4.8) 

where Pj = fulfil \\Pj\\ = 1 (but are otherwise arbitrary), and where pj G (0,1) is a 

null sequence to be specified later. With this choice, we have ||<y9j|| = 1, and one computes from ( |3.23 ) 


that 


where Mj is the 2x2 matrix 


(4.9) 


with the functions 


Mj^ran = J d9 dp hj^rnn{9, p)g^ {pkj{9, p)) X2,p, {9) X2,p, {v) 

(4.10) 

hj,ii{9,p) = hj^22{9,p) = cosh^ (j + ^)Fp{9 - p), 

(4.11) 

hjp2{9,p) = hjpi{9,p) = cosh^ {^)Fp{2j + 9 + p), 

(4.12) 

kj{9,p) = 2sinh(j + ^^)sinh(^^). 

(4.13) 

X is even. One of the eigenvalues of Mj is Xj = Mj^n — Mj^ 12 , 

and we will show 


that Xj —)• —oo, proving that (4.7) holds for a suitable choice of Pj. 

To that end, we establish estimates on hjpi, hj^i 2 and kj for 9,r] G [~PjjPj]i that is, in the region 
where the integrand of (4.10) is nonvanishing. First, continuity of Fp and Fp(0) = 1 imply that 
Fp{9 — p) < c + ^ if Pj is small; note that c > ^ enters here. Estimating cosh^ x < 1 + e^^/4 (x > 0), 
we then have in the relevant range for 9, p, 


hj,u{9,p)< (l + 4e2^e^^^)(c+4). 


Further, the growth condition (4.6) implies for j > 0o/2 + 1, 

hj,i2{9,p) > ccosh(2j + 9 + p) > 


(4.14) 


(4.15) 


Now (4.14) and (4.15) combine to give 

hj,n{9,p) - hj^i2{9,p) < c + ^ (e^^J(c+ 4) - 2ce"^^^) < 2c - c'e^^ < 0 

^ '-V-^ 

->• |-c<0 

with some constant c' > 0 and for large j. Finally, since |0 — r/| < 2pj < 2 in the integrand, 


\kj{9,p)\ < 2{e^e 


< 12e^ Pj. 


(4.16) 


(4.17) 


Now setting specifically pj = ae P with a > 0 still to be specified, we have \kj{9, p)\ < 12a independent 
of j. Noting that 5^(0) = f dtg^{t) > 0, we can achieve with a suitable choice of a that 


1 


g‘^{pkj{9,p)) > -5^(0) > 0. 


(4.18) 


Using (4.16) and (4.18) in the integrand of (4.10), we then obtain 

1 X. /,, ,, N 2 


Xj = Mjpi - Mjp 2 < {2c- c'e‘^^)-g‘^{0 ){\\X2 ,pj\\i) =-g‘^{0){2ce ^ - c'e+^)(p^- ^ ||x 2 ,p,||i) • (4.19) 

_-j /o 

Herep. ||x 2 ,p,||i is actually independent of pj. Hence Xj —>■ —oo as j —)> oo, which concludes the 


proof. 


□ 
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5 Quantum energy inequalities 


We now turn to the existence of quantum energy inequalities, i.e., we want to show that the operator 
is bounded below at one-particle level. As we have seen in Proposition |4.2[ this can be true 
only if the function Fp does not grow too fast. The main goal of the section is the following theorem, 
which establishes a QEI under certain bounds on Fp. 


(5.1) 


Theorem 5.1. Suppose that there exist constants 9 q > 0, Aq > 0, and 0 < c < ^ such that 
|Tp(C)| < c coshReC whenever \Re C\ > do, |IniC| < Aq. 

Further, let g G 5r(]R). Then, there exists Cg > 0 such that 

Vy. G V{R) : > -Cg\M\ (5.2) 

The constant Cg depends on g (and on Fp, hence on P and S) hut not on (p. 


The idea of the proof is as follows. In (3.23), we split the integration region in both 9 and rj into 
the positive and negative half-axis. Setting if {9) := {ip{9),(p{—9))'^ for 0 > 0, we can rewrite our 
expectation value as 


:= {ip,T^^{g^)f) = ^ I d9 I dgg^{gcoshd - gcoshrj)f (9)*M{9,g)if{g), (5.3) 

Jo Jo 


where the matrix M{9, g) is given by 

(cosh? ^JJ^Fp {9 — g) cosh? ^?^Fp {9 + g) 


M{9,g) = 


\cosh? Fp{9 + g) cosh? ^JJ^Fp{9 — g) 


(5.4) 


The eigenvectors of M are Vp = (1,1)^ and v- = (1, —1)^, independent of 9, g, and the corresponding 
eigenvalues are 


h±{9, g) = cosm 


-Fp{9 — g) ± cosm 


-Fp{9 + g). 


2 ^ 2 
Denoting by f± the components of f in the direction of v±, we thus have 

,2 /*c>o /*oo ^ __ 

dg g'^ {fjL cosh 9 — g cosh g) E f±{9)h±{9,g)f±{g). 


^ 2tt 


d9 


(5.5) 


(5.6) 


± 


We will compare to the following related integral expression: 




dg g'^ (/i cosh 9 — g cosh g) E f±{9)k±{9)k±{g)f±{g), 


± 


where 


k±{9) := ^\h±{e,9)\ = ^J\ cosh^ 9 ± Fp{29)\. 


(5.7) 


(5.8) 


Specifically, we will show that > 0 and that \X^p — Yg,\ is bounded in ||¥?|p. We do this in several 
steps; the hypothesis of the theorem is always assumed. 


Lemma 5.2. For any f, we have > 0. 
Proof. Using the identity 


dq 


g^ip-p')= / ^giq + p)9iq + p'), 


we can rewrite the integral as 

,2 


Y — ^ 

^ ~ 47r2 


J dq J dga±{g,q) with a±{g,q) := k±{g)f±(g)g{q + gcoshg). 


But this is clearly nonnegative. 


(5.9) 

(5.10) 

□ 
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For estimating \X^ — y^|, we first need an estimate for the relevant integral kernels, into which 
the growth bound (5.1) for Fp will crucially enter. 


Lemma 5.3. Set 

'■= ^±(p+ i,p- i) - k±(p+^)k±(p - I). (5.11) 

Then, there exists o > 0 such that for all p > Oq + 1 and r G [—1,1], 

K±(p, t) I < ar^ cosh^ p. (5.12) 

Proof. One notes that i±{p, 0) = 0 and that i±{p, r) is symmetric in r. A Taylor expansion of i±{p, r) 
in r around r = 0 then yields that 


\^±{p,r)\ < V sup 
^ kl<i 


d^l 


± 


dr'^ 




(5.13) 


Thus our task is to estimate the derivative. As a first step, we remark that Cauchy’s formula allows 
us to deduce estimates for the derivatives of Fp from (5.1): One finds constants c',c" > 0 such that 


\/ 6 > 0 o + l: 


dFp 


dO 


Now we explicitly compute 

52 


{6) <c'cosh6, 


■^h±{p+ |,p - I) = (^cosh^ pFp{t) ±cosh2 ^Fp 


dSFp 


de^ 


( 0 ) 


< d' cosh 0 . 


(5.14) 


9r2 

d'^Fp 1 

= cosh 2 p-^(r) ± -coshrFp( 2 p). 


(5.15) 


In the region |t| < 1, p > 6*0 + 1, we therefore have due to (5.1), (5.14), 

d 


Qj.2 k±{P + 2 ^ P 2 ^ 


< oi cosh^ p 

with some ai > 0. For the derivative of the second term in ^±, we obtain 


|^fe±(p+ p/c±(p- i) = ^^(p+ |)A:±(p- I) + l;c±(p+ i)^(p- i). 


(5.16) 


(5.17) 


Note here that the radicand in k±{6) - cf. (5.8) - is actually positive for 6 > 6 q due to (5.1), thus the 
function is differentiable, (c < ^ enters here.) For 6 > Oq, we then find \k±{6)\ < 02 cosh0. For its 
second derivative, we have 


= dt 

dO^ 


{e) = ^\l^o^^^e±Fp{2e) 


(4sinh(20)±^^(20))2 ^ cosh( 20 )± 2 ‘|^( 20 ) 


(5.18) 


+ 


(cosh2 e ± Fp(20))3/2 (cosh2 e ± Fp( 20 ))V 2 ' 


In the denominator, we can estimate cosh^ 6 ± Fp{26) > 03 cosh^ 9 (for 9 > 9q, with 03 > 0); again, 
c < ^ enters. In the numerator, we estimate the derivatives of Fp by (5.14). This yields 


d?k 


A 


d 02 


< 04 cosh 9. 


Combining ( |5.16| ), ( |5.17| ) and ( |5.19[ ), we find for p > 0o + 1; kl < Ij 

d'^i± 


9r2 


ipx) 


< as cosh^ p 


with some 05 > 0. Applying this estimate in (5.13) yields the desired result. 


(5.19) 

(5.20) 

□ 
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This estimate on the integral kernel allows us to hnd bounds for ~ T^|- 
Lemma 5.4. For every g G 5]r(]R), there exists Cg > 0 such that —Ty,! < Cg||(/7|p for any ip G 2?(M). 
Proof. Comparing (|5.6p and (5.7) and using Cauchy-Schwarz, we can choose 

|2 


Cn . - 


-E 


poo poo 

/ 

lo Jo 


dg \g'^{pcosh9 — pcos]ig)\ \h±{6,g) — k±{6)k±{g)\" 


(5.21) 


if we can show that these integrals converge. We change to new variables: p = t = 9 — g, 

\d{p,T)/d{9,g)\ = 1. The integral now runs over 0 < p < oo, —2p < t < 2p. However, the integral 
over the region 0 < p < 0o + 1 certainly exists since the integrand is bounded there. Also the integral 
over p > 9 q + 1, 1 <\t\ < 2p is finite, since h±, k± are bounded by a polynomial in coshp, whereas 


|p2(pcosh0 — pcoshp)| = |p2(2psinhpsinh I 


(5.22) 


decays faster then any power of coshp there, as g is of Schwartz class. Thus it remains to prove 
convergence of 

/•r>n /*! 

i2 I ^ / Xi2 


/ dp dr |p2(2psinhpsinh ^)| |^±(p, r)|^ 
' 00+1 • 1-1 


with £± as in Lemma 5.3 We estimate (5.22) more carefully for p>0o + l)|i‘|<l- 

b 


I (p cosh 9 — p cosh g) | < 


cosh^ p + 1 


with some 5 > 0. Using this and the estimates from Lemma 5.3, we have 

roo /•! 


dSlMI ) < / dp 

J 6o-\-l 


r Pf- r < oo. 

_i (r^ cosh^ p + 1)2 Jeo+i cosh p {k^ + 1)2 


dr- 


(We have used the substitution n = r coshp.) This means that Cg in (5.21) is finite. 

Our main result is now a direct consequence of the lemmas above. 

Proof of Theorem\5.1\ Using Lemmas |5.2| and |5.4[ we can compute 


(5.23) 

(5.24) 

(5.25) 

□ 


(p, r‘J"(p2) p)=Y^ + {X^ - Yg,) >Y^- \X^ - Y^\ > -Cg 


□ 


(5.26) 


Note that we have actually shown that T*^®(p^) is a sum of a positive quadratic form and a Hilbert- 
Schmidt operator, since we estimated jXg, — y^| by showing that the integral kernel involved is an 
function. But we will not need this detail for our present purposes. 


6 Uniqueness of the energy density in certain models 

Our results so far indicate that the existence or non-existence of quantum energy inequalities at one- 
particle level is deter mine d by the asymptotic behaviour of the function Fp{9) = P(cos hd)i Tiin (d+i7r). 
Essentially, Theorem 5.1 says that QEIs hold if Fp{9) < ^ coshd, and Proposition |4 . 2| shows that QEIs 


do not hold if Fp(9) > ^ coshd. Let us now see what this means in specific models. 
Let us first consider the Ising model, with the function Emin given by 

Emin(C) = -isinh^. 


( 6 . 1 ) 


If now P = 1, then Fp{9) ~ cosh^'^^ 0 for large 9 (also on a small complex strip), the conditions 
of Theorem |5.1| are fulhlled, and the energy density satisfies a QEI. If, on the ot her hand, P is a 
polynomial of degree d > 1, then Fp{9) ~ cosh'^^^'^^ 9, the conditions of Proposition 4.2 are met, and 
hence no QEI can hold. In other words, in the Ising model, the one-particle energy density is uniquely 
determined by Proposition |3.1| and the additional requirement that a QEI holds. We summarize this 
in the following theorem. 
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Theorem 6.1. In the Ising model (cf. Tafele a one-particle quantum energy inequality holds if, 
and only if, P = 1 in Proposition \3.1\ □ 


We remark that in the case where no QEI holds, Proposition 
unbounded below for any real-valued Schwartz function g ^ 0 - 
restricting to a smaller class of smearing functions. 


4.2 


actually implies that T^^{g 


IS 


we cannot even establish a QEI by 
In the sinh-Gordon model, the situation is more complicated. Here the function Emin is of the form 


Emin(C) = exp Jb(C). 


( 6 . 2 ) 


Since JB(0 + iA) converges to a constant as 0 —>■ ±oo, uniformly on a strip around A = vr, we know that 
|Emin(^+*A)| is bounded the re an d does not converge to 0 at infinity. If now P = 1, then Fp is bounded 
on the strip, thus Theorem 5.1 applies and yields a QEI. If degE > 2, then \Fp{9)\ rsj (cosh(9)^®e^ 


grows faster than ccosh0 for any c, and by Proposition 4.2 no QEI can hold. If however degP = 1, 
then details of the polynomial matter. Since E(l) = 1, such P must be of the form 


P{x) = {1 — i^) lyx with some i/ G M. 

The existence of a QEI now depends on u. Setting := lim 0 _^±oo F m\n {9 + in), we have 

,Fm\n{0 + in) 


Fp{6) = (1 - i')Fmin{9 + in) + iyF ^^„ cosh6>- 


_poo 


(6.3) 


(6.4) 


The hrst summand is negligible against cosh0 for large 9, and the last fraction becomes arbitrarily 
close to 1. Therefore, if |i^| > 1/2E^^, then Proposition 4.2 applies with c = \\v\F^i^ + 2^5- 


No QEI can hold. If, on the other hand, \v\ < 1/2FP^, then similarly Theorem 5.1 applies with 
c = + i ^ follows. In the borderline case \v\ = l/2E^jj, our methods do not 

yield a result. 

Interestingly, the behaviour of the free theory, where Emin = Ij exactly matches the one of the 
sinh-Gordon model, only with ET^ = 1. We summarize: 

Theorem 6.2. In the free and the sinh-Gordon models (cf. Ea5Ze[^, a one-particle quantum energy 
inequality holds if 

1 


P{x) = (1 — v) + vx with 1/ G M, \v\ < 


2 poo ’ 


(6.5) 


where E(^„ = lim 5 i_ 5 .j-oo Emin(^ + ^7r). If P is of this form but with \iy\ > 1/2 E(^„ , or if degP > 2, then 
no quantum energy inequality holds. □ 

Thus, in the sinh-Gordon model and the free field, demanding a QEI strongly restricts the choice 
for the energy density, but it does not hx the energy density uniquely. Even at one-particle level, 
a certain arbitrariness remains. At least in the free field situation, we observe however that the 
“canonical” case E = 1 is distinguished by the fact that it is the only choice which makes the energy 
density nonnegative in one-particle states. 

The generalized Ising and generalized sinh-Gordon models in Table show a similar behaviour, 
but allow for some more variety, since the factors —isinh((C/2) in the denominator of Emin make 
the function decay at real infinity. In consequence, for degE not too large, QEIs will hold due to 


Theorem 5.1; in fact the hypothesis of Proposition 4.1 may be violated and the energy density may 


be positive at one-particle level. For sufficiently large degE, however, Proposition 4.2 will apply and 
hence QEIs will break down. 


7 Numerical results 

In the previous sections, we have discussed whether expectation values of the energy density can be 
negative, and we established under which conditions there are lower bounds to the energy density, i.e., 
T^^{g‘^) > —Cgl at one-particle level. Glearly it would be of interest to compute what the best lower 
bound Cg is for given g, or in other words, what the lowest spectral value of the integral operator in 
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question is. However, all our methods so far are based on estimates, which are in fact quite inexplicit 
at places, and are therefore not suited to compute a useful value for the best lower bound. The 
only feasible way to do so appears to be numerical approximation. In this section, we will describe 
numerical methods to that end. 

To be precise, we consider the operator := PiT^^{g‘^)Pi, where Pi is the projection onto 

the one-particle space. Our task is to compute the lowest value in the spectrum of T^neid'^)- ( So fa r, 
we had defined T^^eid'^) as a symmetric operator on P(]R) only. However, since by Theorem |5.l| it 
is bounded below - we assume that the hypothesis of the theorem is fulfilled -, canonically 

extends to a selfadjoint operator with the same lower bound |RS75[ Thm. X.23], denoted by the same 
symbol, so that has a well-defined spectral decomposition.) 

As a first complication, while we have shown that T^neid'^) i® bounded below under certain con¬ 
ditions, it is certainly unbounded from above due to the strong growth of the integral kernel at large 


arguments, cf. Eq. (3.23). We need to remove this discontinuity before the operator becomes amenable 
to a stable numeric treatment. To that end, we introduce a “rapidity cutoff”: we will restrict the 
support of the wave functions to an interval [—R,R], with sufficiently large P > 0; in other words, 
we consider as an integral operator on L'^{[—R,R\,d9) rather than on L‘^(R,d9). On this 

restricted space, the operator is bounded, in fact it is compact. The cutoff will certainly affect the 
large positive spectral values of the operator; but we claim that it does not significantly change the 
lowest spectral value, which we are interested in here. That this claim is at least self-consistent will 
become clearer below. 

As a next step, we discretize the integral operator. To that end, we use an orthogonal Galerkin 
method with piecewise constant bump functions and application of the midpoint rule; alternatively 
speaking, the Nystrom method associated with the midpoint rule. For details of these approximation 
algorithms see, e.g., |CC831 Ch. 4]; but for the beneht of the reader, let us give a rough description of 
the method here. We divide the interval [—R, P] into N equally sized subintervals of length h = 2R/N, 
with midpoints 9j = —R [j + \)h, where j = 0,..., X — 1. We consider the orthonormal system of 
step functions supported on these intervals, 


ipj{9) = 



if \9j -9\< 
otherwise 


0 <j < N. 


(7.1) 


Then, we will only consider the matrix elements 






2TTh 


d9 


Q _j_ j-j 

dtj cosh^ — Fp{9 — r])g‘^{g,cosh9 — /rcoshry). (7.2) 




For easier evaluation, we will approximate the integral using the midpoint rule, which yields 

2 T 


M, 


jk 


Mjk '■= — hcosh^ 
^ 27r 


Fp{9j — 9 k) g^{g cosh 9j — g cosh 9k). 


(7.3) 


This expression is then straightforward to compute, although evaluating the function Fp will in some 


cases need numerical integration techniques for approximating Jp, cf. (2.5). The eigenvalues and 


eigenvectors of the real symmetric N x N matrix M can now be found with standard numerical 
methods, such as the implicit QL algorithm. While we have not given full error estimates here, we 
can reasonably expect that the lower eigenvectors and eigenvalues of M give an approximation of the 
corresponding properties of T^^^^g"^). The computer code implementing the method, which was used 
to produce the results in the following, is provided with this article [Cod| . along with documentation. 

For the sake of concreteness, we will in the following fix the smearing function 17 to be a Gaussian, 


g{t) = TT 



exp 


8a^) 


with a dimensionless parameter cr > 0 , so that f g{t)‘^ dt = 1 and 


g‘^{p) = exp ( - 


(7.4) 


(7.5) 
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(a) Ising model 



e 


(b) Sinh-Gordon model with coupling B = \. 


Figure 1: Lowest eigenvector of Parameters: N = 500, R = W, a = 0.1, P = 1. 


With this, it is clear that the matrix elements (7.2) scale quadratically with /r, so that we only need 
to consider the “normalized” case fi = 1; in other words, all energy density values will be plotted in 
units of g? in the following. 

We will first discuss the form of the eigenvector for the lowest eigenvalue of M, which turns out to 
be clearly separated from the next higher eigenvalue. Figure shows the eigenvectors (eigenfunctions) 
for the Ising model and for the sinh-Gordon model with maximal coupling (i.e., B = 1), both with 
the “canonical” form of the energy density {P = 1). While the vector certainly differs between these 
models, we can observe in both cases that the eigenfunction decays rapidly for large 9, much before 
the chosen cutoff (here R = 10) is reached. Hence our assumption that the lowest eigenvector is not 
influenced by the (reasonably large) cutoff is self-consistent. 

This observation also yields a useful consistency check for the numerical method: By testing 
whether the components of the eigenvector in direction ^pj are small for j ~ 0 and j ~ — 1 , we can 

verify whether our choice for the cutoff R was appropriate; see the method 

kernels. OneParticleKernel. checkPlausibilityO in the computer code. We remark that a suit¬ 
able choice of R is crucial to the stability of the numeric method: The highest (positive) eigenvalues 
of M are several orders of magnitude larger than the lowest (negative) one, and they grow rapidly 
with increasing R. Hence if R is chosen unduly large, significant roundoff errors will affect the result 
of the QL algorithm when working at a given floating point precision. 

Next, let us compare the lowest eigenvalue of M in various models. Specifically, let us choose 
the sinh-Gordon model, depending on the coupling constant 0 < H < 2, again with P = 1. Figure 
2a shows the lowest eigenvalue in this model as a function of the coupling constant B. Note that B 
and 2 — B correspond to the same scattering function. As expected, when B approaches the values 
0 or 2, the lowest eigenvalue approaches the one of the free theory (namely, 0). The lowest possible 
negative energy density is reached when the coupling is maximal {B = 1), which fits with the picture 
that negative energy density in one-particle states is an effect of self-interaction in the quantum field 
theory. 

We also investigate the corresponding generalised Ising model with one coupling constant, i.e., 
with scattering function 


5(c)=- ™;;):‘""n"(^ -WGo,.o„(c). o<b<2. 

smh C + * sm Ptt /2 


(7.6) 


Note that S{() —>• S'ising(C) = — 1 as P —)■ 0 . In this case, rather than P = 1 we choose P{x) = (l-|-x)/ 2 : 
this guarantees that also Fp converges to the canonical Ising model expression as P —)• 0. Both 


Proposition |4.1| and Theorem 5.1 are applicable with this choice of P. As Fig. |2b| shows, the lowest 
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(a) Sinh-Gordon model, coupling 0 < -B < 2 (b) Generalised Ising model, coupling 0 < B < 2 

Figure 2: Lowest eigenvalue of depending on the coupling constant. Parameters: N = 500, 

R = 7, a = 0.1; (a) P = 1, (b) P{x) = (1 + x)/2. 


possible eigenvalue is still reached at maximal coupling B = 1, whereas for B —)• 0 we obtain the Ising 
value of Ri —0.116. 




(a) Several real coupling constants with B = 1 


(b) Several complex coupling constants Bk = l + k-O.li, 
1 < k < n, with complex conjugate factors added. 


Figure 3: Lowest eigenvalue of Fone(5'^) the generalized sinh-Gordon model. The horizontal axis 
shows the number of coupling constants n. Parameters: N = 500, R = 7, a = 0.1; see Eq. (7.7) for 
P. Note the logarithmic scale on the vertical axis. 


To investigate other types of interaction, we also show the lowest eigenvalue in the generalized 
sinh-Gordon model with several factors in the scattering matrix, i.e., with several coupling constants 
Bj (cf. Table [^. We do this once for n coupling constants that are equal and real {Bj = 1 for all 
j), which is shown in Fig. 3a, and once for pairs of complex-conjugate factors with B = 1 ± A: • O.li, 
k = 1,... ,n see Fig. |3b[ We choose the polynomial P such that the free-held expression Fp = 1 is 
reached when all coupling constants are set to 0; namely 



Ln/2J 


n 


(n real coupling constants) and 
(n complex conjugate pairs). 


(7.7) 
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In both cases, when the number n grows, the lowest eigenvalue of rapidly grows towards 

negative values. This is consistent with our computations in Sec. since in particular the derivatives 
of Fp grow rapidly with n. Thus again, a stronger interaction leads to more negative energy density. 



Figure 4: Lowest eigenvalue in the sinh-Gordon model {B = 1) depending on the smearing width 
parameter a. Parameters: N = 500, R = 10, P = 1. Note the double logarithmic scale 


So far, we have always used the same “smearing function” in our analysis. We now vary the 
width a of the Gaussian function g, see Eq. and observe the change in the lowest eigenvalue of 

the smeared energy density, as shown in Fig. Since the energy density is expected to be a quantum 
field of scaling dimension 2 in the high energy limit, one conjectures that for small values of a, the 
lowest eigenvalue of Tone(5^) roughly scales with In fact, this is confirmed by our numerical 

results; a corresponding linear asymptote with slope 2 has been added to the double-logarithmic plot 
Fig-i to demonstrate this. 



(a) Free field (b) sinh-Gordon model, B = 1 

Figure 5: Lowest eigenvalue of T^^g{g‘^) with a “non-canonical” energy density with P{x) = — 

Parameters: N = 300, i? = 8, a = O.I. The vertical lines indicate the theoretical limit for existence 
of QEIs, \iy\ = l/2F^j^. Note that in case (b), we have 1.267. 

Finally, in the free model and the sinh-Gordon model, we want to investigate the “noncanonical” 
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form of the energy density as described in Sec. where the polynomial P is not constant but linear, 


see Eq. (6.5). In this case, we know from Theorem 6.2 that is bounded below if |i^| < 1/2E‘ 


7100 

min’ 


but unbounded if > 1/2F^^^. Our numerical results in Fig. [^confirm this: At the critical value 
\iy\ = 1/2F^^^, indicated by the vertical lines, the lowest eigenvalue of M becomes several orders of 
magnitude larger than for smaller values of |z/|. In fact, for \i'\ > 1/2E^^, the form of the eigenvector 
changes drastically from what is known from Fig. and the eigenvalue is highly dependent on the 
value of the cutoff R. (Details of this are not shown in the graph.) 


8 Conclusions and outlook 


In this article, we have investigated the energy density in certain integrable models at one-particle 
level. Starting from the general properties of a local operator, as derived in |BC15j . and the generic 
properties of a stress-energy tensor, we were able to determine in one-particle matrix elements up 
to a certain polynomial factor, P(cosh(0 — r/)). While one-particle states with negative energy density 
quite generally exist, the choice of energy density is further restricted if we demand that a quantum 
energy inequality holds. 

In fact, in the Ising model, the existence of one-particle state-independent QEIs fixes the form 
of the energy density in one-particle states uniquely. In a different class of models, including the 
sinh-Gordon model and the free Bose field, QEIs reduce the ambiguity in the energy density, but they 


still leave us with one parameter to be chosen (see Theorem 6.2). This choice affects the form factor 


F 2 of the stress-energy tensor, and is different from the ambiguity in Fi discussed in |MS94j . which 
does not contribute to one-particle expectation values. 

It is interesting to note how this choice in F 2 works out in the free field case {S = 1). The usual 
“canonical” stress-energy tensor. 


^ canonical 


1 //2 

= :d°‘(j) (j)-. — :d'^(j) d<y(f>: + , 


( 8 . 1 ) 


corresponds to our choice P{x) = 1 at one-particle level. But other choices of P can be linked with local 
fields on Fock space as well. In particular, one possible choice of stress-energy tensor corresponding 
to P{x) = X is 


pa/3 

^ alternative 


= -'48^ 


: -h -.8'^(f)8^4 - . 


( 8 . 2 ) 


One can easily show that this field does in fact satisfy the standard properties of a stress-energy tensor 
on the full Fock space, including covariance and the continuity equation. A permissible energy density 
that fulfils QEIs at one-particle level would then be given by 


= 0 — lA -I- 

1 / v ^ canonical “r j- alternative 


with 


1 

< 2 - 


(8.3) 


We have not verified at this point whether the QEI for this energy density holds on all Fock space, 
however. 

Throughout this article, we have considered QEIs only on the level of one-particle states. A 
natural next step would be to see whether, in the same class of integrable quantum field theories, 
this inequality holds for all (suitably regular) states, independent of the particle number. We indeed 
expect that a state-independent QEI will still hold in this case; numerical evidence in two-particle 
states in the sinh-Gordon model suggests this. But establishing such a QEI rigorously will require 
substantially more technical effort. 


Specifically, let us reconsider the expansion (2.2). While in the one-particle case, only the kernel 
E 2 [A] contributes, in the general case one needs to consider all Efc[A]. In order to establish QEIs, we 
would need to compute operator norm bounds on these (or related) integral operators. However, while 
E 2 [A] is analytic, the higher-order Efc[A] will possess (regularized) first-order poles in the integration 
variable, i.e., they are “singular integral kernels”. Methods for finding operator bounds on these 
singular integral operators exist |Ghr90] . but they are much more intricate than those we used in 
Sec. where we essentially estimated the kernel itself in norm. 
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Independent of that, it may also be worthwhile to generalize our results to a larger class of integrable 
models: While in the present paper we have considered one species of scalar Bosons, one would also 
like to consider models with several particle species, possibly of different masses and of higher spin, 
models with inner symmetries (gauge symmetries) such as the non-linear 0{N) a models |BFK13j . 
and models with bound states such as the sine-Gordon model |BFKZ99] . The construction of models 
with several particle species has been investigated in the algebraic approach to integrable systems 
iLsn], which can provide a starting point for our investigation. For models with bounds states, a 
rigorous treatment has only partially been achieved (see |CT15) ). Extending our analysis to these 
models, one might test whether the QEI (1.1) persists for more general types of self-interaction, and 
one would be able to study the behaviour of the lower bounds on a wider range of parameters. 
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